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SUMMARY 


This  report  covers  two  aspects  of  the  short  wave  regime  of 
planetary  radiation.  First,  extensive  measurements  of  the  reflection 
properties  of  a  number  of  surfaces  have  been  made,  both  the  directional 
reflectance  and  the  degree  of  polarization  of  the  reflected  radiation  having 
been  determined.  It  is  found  that  most  surfaces  show  pronounced  "limb- 
brightening"  for  large  angles  of  incidence,  while  for  small  angles  of  inci¬ 
dence  the  reflectance  does  not  vary  greatly  with  angle  at  which  the  surface 
is  viewed.  The  degree  of  polarization,  shows  a  characteristic  pattern  for 
most  surfaces,  a  maximum  being  located  in  the  principal  plane  at  100°  to 
120°  from  the  antisource  direction  and  a  small  region  of  negative  polariza¬ 
tion  being  observed  in  the  vicinity  of  the  antisource  direction.  For  most 
surfaces  the  reflectance  increases  and  the  maximum  degree  of  polorization 
decreases  with  increasing  wavelength. 

The  second  part  of  the  research  has  been  the  introduction  of  the  re¬ 
flectance  values  into  the  theory  of  radiative  transfer  to  determine  the  char¬ 
acteristics  of  the  radiation  field  at  high  altitudes  above  the  surface.  It  is 
seen  that  fc  r  the  cases  of  clear  and  slightly  turbid  atmospheres  the  surface 
properties  have  a  strong  influence  on  the  contrasts  which  would  be  observed 
from  a  high  altitude  vantage  point.  A  possibility  of  using  polarization  char¬ 
acteristics  to  enhance  contrasts  is  demonstrated. 
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Abstract 


This  report  covers  two  aspects  of  the  short  wave  regime  of  planetary  radia¬ 
tion.  First,  extensive  measurements  of  the  reflection  properties  of  a  number  of 
surfaces  have  been  made,  both  the  directional  reflectance  and  the  degree  of 
polarization  of  the  reflected  radiation  having  been  determined.  It  is  found  that 
most  surfaces  show  pronounced  "limb-brightening"  for  large  angles  of  incidence, 
while  for  small  angles  of  incidence  the  reflectance  does  not  vary  greatly  with 
angle  at  which  the  surface  is  viewed.  The  degree  of  polarization  shows  a 

characteristic  pattern  for  most  surfaces,  a  maximum  being  located  in  the 

o  o 

principal  plane  at  100  to  120  from  the  antisource  direction  and  a  small  region 
of  negative  polarization  being  observed  in  the  vicinity  of  the  antisource  direction. 
For  most  surfaces  the  reflectance  increases  and  the  maximum  degree  of  polariza¬ 
tion  decreases  with  increasing  wavelength. 

The  second  part  of  the  research  has  been  the  introduction  of  the  reflectance 
values  into  the  theory  of  radiative  transfer  to  determine  the  characteristics  of 
the  radiation  field  at  high  altitudes  above  the  surface.  It  is  seen  that  for  the 
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strong  influence  on  the  contrasts  which  would  be  observed  from  a  high  altitude 
vantage  point.  A  possibility  of  using  polarization  characteristics  to  enhance 


contrasts  is  demonstrated. 
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I.  Introduction 


It  is  well  known  that  the  radiation  which  is  directed  to  space  from  the  upper 
limits  of  the  atmosphere,  as  well  as  that  which  is  proceeding  upward  inside  the 
atmosphere,  is  composed  partially  of  radiation  which  has  been  reflected  from  the 
planetary  surface.  One  has  only  to  look  at  television  pictures  taken  from  the 
TIROS  or  Nimbus  meteorological  satellites,  at  photographs  obtained  by  Astro¬ 
nauts  while  in  orbit  around  the  Earth,  at  photogiaphs  taken  on  an  aerial  mapping 
or  reconnaissance  mission,  or  indeed  at  the  sunlit  scene  presented  at  the  win¬ 
dow  of  a  high  altitude  aircraft,  to  see  just  how  important  surface  reflection  is 
in  determining  the  characteristics  of  the  upwelling  radiation.  The  intensity  and 
polarization  of  the  radiation  observed  in  such  circumstances  are  obviously 
functions  of  the  angle  at  which  the  sunlight  strikes  the  surface,  the  angle  at 
which  the  surface  is  viewed,  the  wavelength  of  the  radiation,  the  physical  proper¬ 
ties  of  the  surface  itself,  and  the  transmission  and  scattering  properties  of  the 
intervening  atmosphere.  The  present  research  work  has  been  an  investigation 
of  the  effects  introduced  by  these  various  parameters  individually  and  of  the 
characteristics  of  the  total  radiation  which  proceeds  outward  through  the  upper 
reaches  of  the  atmosphere. 

Measurements  by  different  authors  have  already  indicated  some  of  the 
reflective  properties  of  various  types  of  surfaces.  The  results  obtained  by 
Miller  (1955)  show  that  a  new  snow  surface  may  reflect  more  than  80  percent 
of  the  visible  radiation  which  is  incident  on  it.  The  extensive  measurements 
by  Krinov  (1947)  show  total  reflectance  values  of  soil  and  rocks  varying  from 
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as  little  as  0,  02  for  black  candy  loam  to  0.  35  -  0,  75  for  highly  reflecting  clay, 
limestone,  and  shale,  Bauer  and  Dutton  (1962)  measured  wavelength-integrated 
total  hemispheric  reflectance  of  farmlands  and  wooded  hills  of  Wisconsin  by 
radiometers  mounted  on  aircraft,  obtaining  values  of  0.  10  to  0.  20  in  the  absence 
of  snow  cover.  Snow  increased  <he  average  reflectance  to  0,  50  -  0.  80  for  other¬ 
wise  similar  conditions,  Except  for  the  effects  of  snow,  the  measurements 
showed  little  change  of  total  reflectance  with  season,  although  it  is  to  be  expected 
that  the  wavelength  distribution  of  the  reflected  radiation  would  vary  with  seasonal 
changes  cf  vegetation. 

The  variation  of  reflectance  with  angle  at  which  radiation  is  incident  on  the 
surface  and  angle  at  which  the  surface  is  viewed  has  received  much  less  attention 
than  has  the  total  hemispherical  reflectance.  The  effect  of  angle  of  incidence  has 
been  investigated,  both  theoretically  and  experimentally,  more  for  water  sur¬ 
faces  than  for  other  types  of  surfaces.  Of  particular  note  here  are  the  works  of 
Cox  and  Munk  (1955)  in  determining  the  effects  of  surface  roughness  on  the 
reflectance  properties  of  the  sea,  and  of  Anderson  (1952)  in  measuring  the  energy 
reflected  from  an  inland  lake.  The  theoretically  predicted  strong  dependence  of 
reflectance  of  water  on  incident  angle  is  borne  out  in  the  measurements,  although 
roughness  of  the  sea  surface  is  influential  in  determining  both  the  total  reflec¬ 
tance  and  the  intensity  distribution  of  the  reflected  light.  For  natural  land 
surfaces,  the  effect  of  incident  angle  is  considered  by  Kondratev  and  Manolova 
(1955)  in  relation  to  the  radiation  balance  of  slopes,  Ashburn  and  Weldon  (1956) 
and  Coulson  (1956)  showed  a  considerable  dependence  of  the  reflection  properties 
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of  semi-desert  surfaces  on  solar  zenith  angle,  and  a  number  of  investigators, 
Middleton  and  Mungall  (1952),  and  Christie  (1953),  for  example,  have  shown 
a  strong  functional  dependence  of  snow  reflection  on  angle  of  incidence. 

Measurements  of  the  directional  reflectance  properties  of  terrestrial 
surfaces,  as  evidenced  by  the  angular  intensity  distribution  of  the  reflected 
light,  have  been  oriented  mainly  toward  interpreting  the  reflectance  properties 
of  the  Moon.  Orlova  (1956)  has  made  reflectance  measurements  of  a  number  of 
surfaces  in  the  laboratory  and  classified  the  "indicatrices"  into  orthotropic, 
reflecting,  rough,  and  mixed  categories,  depending  on  the  number  and  position 
of  reflection  maxima.  Additional  measurements  on  very  complex  terrestrial 
materials,  similarly  directed  toward  the  interpretation  of  photometric  properties 
of  the  Moon,  have  been  made  by  vanDiggelen  (1959),  Hapke  and  vanHorn  (1963), 
and  others.  Most  of  these  measurements  have  been  confined  to  the  principal 
plane.  Krinov's  (1947)  measurements  were  made  mainly  in  the  direction  of  the 
nadir,  but  for  a  few  cases  the  instrument  was  directed  at  a  45°  nadir  angle  and 
90°  azimuth  angle  with  respect  to  the  azimuth  of  the  sun. 

Measurements  on  the  polarizing  properties  of  materials  extend  back  to 
the  pioneering  work  of  Arago  and  to  the  relatively  extensive  set  of  measure¬ 
ments  performed  by  Brewster  (1865)  on  painted  surfaces,  paper,  cloth,  snow, 
and  white  powders  of  various  kinds.  Lyot  and  Dollfus  (1949),  by  the  use  of  a 
polarimeter  invented  by  Lyot,  determined  the  degree  of  polarization  of  light 
reflected  by  the  Moon.  The  polarizing  properties  of  a  number  of  types  of 
natural  and  artificially-produced  materials  have  been  measured  by  Dollfus 
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(1961,  1957)  vising  the  technique  of  Lyot.  Dollfus  applied  his  polarization 
results  to  interpreting  the  composition  of  the  lunar  and  Martian  surfaces  and 
to  a  deduction  of  an  upper  limit  for  the  density  of  the  atmosphere  of  the  Moon 
and  in  determining  the  total  mass  of  the  atmosphere  of  Mars. 

Unfortunately  none  of  these  investigations  has  been  adequate  for  present 
requirements  in  equipment  development,  more  complete  data  on  all  of  the 
different  parameters  being  required.  The  general  approach  taken  in  this  in¬ 
vestigation  was  designed  to  provide  results  of  maximum  utility  in  the  practical 
problem  of  object  detection  by  an  imaging  device,  be  the  device  the  human  eye 
or  a  photographic  or  television  camera.  In  particular,  the  possibilities  of  using 
the  polarization  characteristics  of  the  radiation  to  enhance  contrasts  in  the 
viewed  field  are  of  major  interest.  The  investigation  has  encompassed  measure¬ 
ments  of  the  directional  reflectance  and  degree  of  polarization  of  the  reflected 
radiation  for  ten  different  types  of  surfaces,  the  instrumental  requirements 
for  making  those  measurements,  the  theory  of  radiative  transfer  by  which  the 
measurements  can  be  quantitatively  interpreted  in  the  context  of  the  problem, 
and  the  computational  results  obtained.  Each  of  these  aspects  will  be  discussed 
below. 

XI.  Instrumentation 

The  instrument  which  has  been  used  for  the  measurements  is  shown 
schematically  in  Fig.  1.  The  material  sample  is  contained  in  a  2  feet  square 
sample  tray  S  which  is  illuminated  by  a  light  source  L  from  a  zenith  angle 


4 


Fig.  I  -  Schematic  representation  of 
apparatus  used  in  reflection  measurement 
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Fig.  2  -  Schematic  diagram  of 

nVintnpl  prf  ri  r*  roflorfnmofor 

r  —  —  —  —  —  — '-'--v-  .1  “ 
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The  surface  is  viewed  from  a  zenith  angle  0  and  azimuth  by  the  photo¬ 
electric  reflectometer  R. 

The  reflectometer  itself  is  shown  schematically  in  Fig.  2.  The  light 
L.  enters  the  collimator  tube  C,  is  chopped  by  a  chopper  rotated  by  motor  M, 
and  passes  successively  through  a  rotating  analyzer  A  and  optical  filter  F,  and 
finally  activates  a  radiation  detector  D.  The  detector  for  radiation  of  wave¬ 
lengths  ^  s  796oX  has  been  a  photomultiplier  tube  (RCA),  and  for  X  =  1.025(1 
a  silicon  photodiode  has  been  used.  The  response  of  the  detector  is  amplified 
and  recorded  on  a  standard  strip-chart  recorder.  The  acceptance  field  for 
the  instrument  is  limited  to  a  2.  5°  half-angle  cone  by  the  diaphragms  within 
the  collimator  tube.  In  order  to  minimize  possible  polarization  sources  within 
the  instrument,  the  optical  components  have  been  restricted  to  the  bare  essen¬ 
tials,  no  lenses  or  mirrors  being  in  the  optical  train.  The  analyzer  which  has 
been  used  for  most  of  the  measurements  is  a  sheet  polarizer  (Polaroid  HN-22) 
which  gives  better  than  99.  9%  polarization  throughout  the  spectral  range  of 
X  ^  7960A,  A  Glan-Thompson  prism  was  used  at  X  =  l.  025  M . 

The  measurements  have  been  made  in  four  spectral  regions  defined  by 
interference  filters  and  centered  at  wavelengths  4920A,  6430A,  79bOA,  and 
1.  025  (i  .  Transmission  curves  for  the  filters  are  shown  in  Fig.  3.  The  first 


three  of  those  particular  wavelengths  ccufs^ouu,  rc&^cv-tiv^ly,  tu  riorxnai 


optical  thicknesses  of  0.  15,  0.05,  and  0.02  for  the  Earth's  molecular  atmosphere. 


The  wavelength  of  1.  025^  was  chosen  so  as  to  minimize  atmospheric  effects. 


both  of  scattering  by  the  atmospheric  components  and  absorption  by  atmospheric 
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water  vapor. 


A  schematic  diagram  of  the  components  ol  the  reflectometer ,  for  the 
case  in  which  the  photomultiplier  tube  was  used  as  the  detector,  is  shown  as 
Fig.  4.  The  basic  configuration  is  a  standard  one  for  a  photomultiplier  tube. 

The  most  stringent  requirement  for  the  present  application  was  that  of  obtaining 
the  required  sensitivity.  A  low  reflectance  of  some  of  the  surfaces  combined 
with  a  narrow  bandpass  filter  resulted  in  a  very  small  energy  flux  incident  on 
the  detector.  The  problem  was  particularly  severe  at  796oX,  for  which  case 
the  S-20  type  cathode  surface  has  a  low  sensitivity.  A  tube  with  an  S-l  response 
was  not  available  in  time  for  use  on  this  short  period  contract. 

These  difficulties  required  a  high-gain  amplifier  to  be  installed  in  the 
system,  a  fact  which  introduced  a  severe  noise  pickup  problem.  It  was  necessary 
to  entirely  repackage  the  amplifier  in  order  to  bring  the  noise  down  to  a  tolerable 
level.  The  repackaging  task  required  a  considerable  expenditure  of  both  time 
and  effort  and  caused  the  measurement  program  to  fall  behind  in  the  schedule 
which  had  been  planned. 

A  concentrated  effort  was  made  to  bring  the  measurements  back  to  schedule, 
a  goal  which  was  only  partially  achieved.  The  net  result  of  the  whole  instrumentation 
problem  was  to  eliminate  the  possibility,  within  the  time  and  funds  available  on  the 
contract,  of  making  measurements  in  the  ultraviolet  spectral  region  as  had  been 
anticipated. 

The  circuit  diagram  of  the  final  electronic  configuration  used  with  the 
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PHOTOMULTIPLIE] 


Schematic  diagram  of  the 
electronic  components  of  the 
reflectom  eter. 


photomultiplier  tube  is  included  as  Fig.  5.  In  making  the  measurements,  a 
procedure  was  established  by  which  an  instrument  calibration  was  performed, 
by  means  of  a  series  of  calibrated  neutral  density  filters,  after  the  instrument 
had  stabilized  at  the  beginning  of  each  day's  measurements.  This  helped  to 
assure  the  validity  of  the  measurements  and  it  turned  out  to  be  very  useful  in 
final  interpretation  of  the  data.  The  system,  response  was  found  to  change  at 
two  different  times  in  the  measurement  program.  At  the  beginning  of  the 
measurements,  the  response  was  linear  with  radiation  intensity  throughout 
the  full  scale  deflection  but  there  was  a  constant  scale  offset  correction  which 
had  to  be  applied.  This  calibration  curve,  which  is  plotted  as  Curve  I  in 
Fig.  6,  was  valid  until  September  15,  at  which  time  the  major  modification 
to  the  amplifier  was  made. 

The  amplifier  modification  introduced  a  nonlinearity  into  the  instrument 
response,  as  shown  by  the  calibration  Curve  II  of  Fig.  6.  This  nonlinearity 
was  taken  into  account  during  the  data  reduction  process  by  fitting  the  curve  by 
a  third  degree  polynomial  and  programming  the  computer  to  automatically  cor¬ 
rect  each  measurement  value  appropriately.  This  added  only  a  negligible  in¬ 
crease  in  the  computer  time  required.  Curve  II  was  applicable  for  the  period 
September  25  to  October  25,  inclusive. 

At  the  start  of  measurement  s  on  October  26  it  was  found  that  a  slight  shift 
of  the  calibration  for  low  light  levels  had  occurred  between  the  calibrations  of 
October  25  and  October  26,  while  the  original  curve  was  still  valid  at  deflections 
greater  than  about  fifty  scale  divisions.  The  new  response,  shown  by  Curve  III 
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in  Fig.  6,  changed  very  little  up  until  the  end  of  the  measurements  with  the 
photomultiplier,  and  it  was  fitted  with  another  third  degree  polynomial  for 
data  reduction.  The  reason  for  the  relatively  sudden  shift  has  still  not  been 
determined. 

The  final  calibration  curve  of  the  equipment,  Curve  IV  of  Fig.  6,  is  that 
applicable  for  the  silicon  detector  and  its  attendant  electronics,  A  slight  non¬ 
linearity  was  again  taken  account  of  by  a  polynomial  approximation  in  the  data 
reduction  program. 

In  order  to  extend  the  spectral  range  of  the  measurements  to  beyond  one 
micron,  it  was  necessary  to  replace  the  photomultiplier  with  a  silicon  photo¬ 
diode  detector  and  to  install  a  Gian- Thompson  awiyzer  in  place  of  the  Polaroid 
disk.  The  photodiode  was  supplied  by  the  Navy  and  the  Gian- Thompson  prism 
was  available  from  another  project. 

Sensitivity  problems  again  were  encountered  and  were  finally  overcome  by 
increasing  the  incident  radiation  by  the  use  of  two  infrared  heat  lamps  as  the 
source.  This  had  an  undesirable  effect  of  increasing  the  dimensions  of  the 
source  itself,  but  it  was  the  only  alternative  available  within  the  scope  of  the 
project.  Increasing  the  area  of  the  source  tends  co  smooth  out  the  angular 
dependence  in  the  measurements. 

Thp  #>1  Artrnnir  ri  rrnit  dii^rEITl  r'£  th°  arnnUfi  gij  US6d  V/it.h  thfi  silicon 

photodiode  detector  is  shown  as  Fig.  7. 
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Fig.  7  -  Infrared  (1.  025  M ) 


Detector  and  AC  Amplifier  Circuit. 
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III.  Measurements 


A.  Theory  of  measurements 

A  radiation  field  in  an  arbitrary  state  of  polarization  can  be  represented 
in  terms  of  four  quantities,  which  are  usually  called  the  Stokes  parameters  after 
the  English  physicist  of  the  19th  century,  Sir  George  Stokes.  If  we  denote  two 
orthogonal  directions  in  the  plane  normal  to  the  direction  of  propagation  of  the 
radiation  by  i  ai,d  j  ,  then  the  intensities  J.  and  L  are,  respectively,  the 
squares  of  the  amplitudes  of  the  electric  vector  along  the  i  and  j  directions. 

The  total  intensity  of  the  radiation,  which  is  the  first  of  the  Stokes  parameters, 
is  given  by 


1=1.  +  I. 

i  J 


(1) 


The  other  three  of  the  Stokes  parameters  are  defined  in  terms  of  the  orthogonal 
components  as  follows: 


U  =  2(Ij  Ij  )1/2cos  6  (2) 

V  =  2(Ii  Ij  f/^in  6 


6  is  the  phase  angle  between  the  vibrations  in  the  i  and  j  directions.  If  0  is 
the  angle  between  the  direction  i  and  the  major  axis  of  the  ellipse  described 
by  the  vibration  of  the  electric  vector  of  the  radiation,  then 

U  =  Q  tan  2  #.  (3) 

The  parameter  V  is  expressed  in  terms  of  the  major  and  minor  axes  of  the 
ellipse,  denoted  by  a  and  b  respectively,  by  the  relations 
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V  =  I  sin  2  0 

B  b  <4> 

P  =  arctan  {  —  ) 

The  degree  of  polarization  is 

.  2  .  2  .  2.1/2 

P  =  (Q  *  f  V  )  (5j 

I 

In  the  present  measurements  it  is  assumed  that  the  only  type  of  polarization 
present  is  linear,  in  which  case 

V  *  o  (6) 

Then  the  only  quantities  necessary  to  characterize  this  partially  linearly 
polarized  radiation  are  I-  ,  Ij  ,  and  ^  .  In  the  reflectometer  being  used,  the 
output  from  the  detector  is  monitored  continuously  as  the  analyzer  is  rotated 
around  the  optical  axis  of  the  instrument.  This  means  that  the  maximum  signal 
is  obtained  when  the  plane  of  transmission  of  the  analyzer  is  parallel  to  the 
direction  i,  in  which  case  4>=  0.  The  minimum  is  obtained  after  another  90° 
rotation  of  the  analyzer,  or  at  0  =  /2.  In  both  cases  U  =  0,  and  the  intensity 

and  degree  of  polarization  of  the  measured  radiation  are  given  by 


1=1.  +  I. 

1  J 

=  I  + 

max 

*min 

r  h  -h 

*max  “ 

*min 

~k 

1  ..  + 

n  i  ,ch  x 

^min 

Hi  riorti nna  1 

*  V*  VU 

Directional  reflectance  of  the  surface  is  obtained  in  the  following  manner. 
It  is  known  that  a  surface  composed  of  magnesium  oxide  represents,  to  a  good 
approximation,  a  perfectly  diffuse  surface  with  a  reflectance  of  1.0  for  visible 
light.  Such  a  surface  is  used  as  a  standard  reflector.  If  F  units  of  energy  per 
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Unit  time  and  unit  frequency  interval  are  incident  on  a  unit  area  of  the  standard 
surface,  the  intensity  IR  of  energy  reflected  into  an  incremental  solid  angle  d  w 
is  independent  of  direction  and  satisfies  the  relation 

lMd“=R6F  =  F,  (8) 


where  p  =  cos  0,  0  being  the  nadir  angle  at  which  the  surface  is  viewed.  Rg 
represents  the  total  reflectance  of  the  standard  surface,  which  is  taken  as  unity 
for  the  magnesium  oxide,  and  the  integration  is  taken  over  the  entire  hemisphere. 
A  directional  reflection  coefficient  also  has  significance,  and  can  be  defined 


by  the  relation 

p  =  ii 

s  F 


(9) 


Obviously  P  is  a  constant  for  a  perfect  standard  (Lambert)  surface.  By  Eq. 
s 

(8)  and  Eq.  (9),  integration  over  a  hemisphere  yields,  for  the  standard  surface, 

(10) 


J  M  d  uJ  =  rs  =  1 

O 


and 


v 


(11) 


Let  I(P,<P)  denote  the  intensity  of  radiation  reflected  from  a  sample  surface 
into  a  direction  given  by  the  nadir  angle  8  and  azimuth  <P  for  the  condition  in 
which  the  sample  is  likewise  illuminated  by  F  units  of  energy  per  unit  time  and 
unit  frequency  interval.  We  can  now  define  a  direction-dependent  reflection 
coefficient  P(K<P)  for  the  surface  by  the  relation 

I  (P,<P)  -  P  )  F.  (12) 
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These  quantities  must  satisfy  the  relations 


I(M,<P)^dto  =  rf 

(13) 

Jft  P(»  ,<P  )*  d  w  =  R  , 

(14) 

where  R  is  the  ratio  of  the  total  reflected  energy  to  the  total  incident  energy. 

The  procedure  followed  for  the  actual  measurements  is  to  illuminate  the 
magnesium  oxide  standard  surface  from  a  direction  given  by  the  zenith  angle 
6Q  £  arccos  HQ  and  azimuth  angle  <P0  =  0  and  set  the  instrument  to  view  the 
standard  surface  from  a  direction  (6,<P).  If  the  standard  were  a  perfect  Lambert 
surface,  the  measured  intensity  would  be  independent  of  0  and  <P .  However,  in 
order  to  assure  uniformity  in  the  measurements,  the  standard  surface  is  always 
'  viewed  from  the  direction  of  the  surface-normal,  in  which  case  6  =  0  and 
is  indeterminate.  Thus,  an  instrument  response  Dg  for  the  standard  surface 
is  obtained.  The  intensity  Is  of  the  reflected  light  is  then 

I5  =  *  Ds  ,  05) 

where  Ag  is  the  instrument  transfer  function  at  intensity  Is. 

The  standard  surface  is  then  replaced  by  the  sample  for  which  the  reflec¬ 
tion  characteristics  are  desired.  The  sample  is  viewed  from  any  desired 
direction,  given  by  (6 and  a  second  instrument  response  D(S,<f$  is  obtained. 
The  intensity  1(6, ip)  is 

I  (6 ,ip)  =  A  D  (6,  <p  )  (16) 

where  A  is  the  instrument  transfer  function  at  this  new  intensity. 

The  instrument  transfer  function  may  or  may  not  be  independent  of  intensity, 
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the  only  requirement  being  that  it  be  a  known  function.  The  most  convenient 
situation  is  that  A  is  a  constant  over  the  range  of  intensities  of  interest,  or 
A  ~  A.  This  ideal  situation  existed  during  the  first  part  of  the  present 
measurements,  as  was  seen  by  the  linear  behavior  of  calibration  Curve  I 
of  Fig.  6  above.  Unfortunately  the  instrument  response  has  been  non-linear 
since  the  modifications  introduced  after  September  15.  A  method  of  introducing 
the  corrections  by  computer  was  devised,  however,  and  the  non-linear  behavior 
has  caused  no  significant  problem  in  the  measurements. 

On  taking  the  ratio  of  Eq.  (16)  to  Eq.  (15)  and  introducing  Eq.  (9)  and 
Eq.  (12)  we  obtain  - 

I  (6  ,<^)  =  D  (B  ,(p)  =  p  (6  ,  <p )  (17) 

i,  (0,-)  D  (0,-)  P TO,"-)” 

s  s  s 

The  assumption  has  been  made  that  the  incident  flux  F  is  constant. 

Finally,  by  introducing  the  known  value  of  P  given  by  Eq.  (11),  we 

s 

obtain  the  final  expression  for  the  directional  reflectance  of  the  sample  surface 
in  the  (0  ,<P)  direction  as 

P(6,<p)  =  1  D  (6  ,<p)  ,  (18) 

V  D  (0. -5 

s 

B.  Mode  of  Operation 

The  mode  of  operation  in  the  measurements  was  typically  as  follows: 

The  sample  was  placed  in  the  sample  tray  in  such  a  manner  that  the 
surface  was  level  with  the  sides  of  the  tray.  For  soil  and  sand  samples  the 
surface  was  carefully  levelled  off  by  means  of  a  straight  edge.  Every  effort 
was  made  to  avoid  directional  effects  being  introduced  during  the  levelling 
process.  For  plant  materials  the  effective  top  of  the  plants,  as  judged  by 


eye,  was  made  to  coincide  with  the  level  of  the  sides  of  the  sample  tray. 

The  surfaces  ox  solid  materials  such  as  concrete  and  blacktop  were  easily 
adjusted  to  the  correct  height,  and  the  sample  of  canvas  was  tightly  stretched 
over  the  top  of  the  sample  tray. 

Once  the  sample  was  properly  prepared,  it  was  placed  on  the  stage  of 

the  reflectometer  mount.  The  stage  is  constructed  so  that  it  rotates  in 

azimuth  as  the  reflectometer  itself  is  rotated.  The  light  source,  a  tungsten 

filament  bulb,  was  then  positioned  for  a  given  angle  of  incidence.  In  these 

o  o 

measurements  three  angles  of  incidence  were  used,  namely,  0  ,  53.1  ,  and 

„  o 

78.  5  . 

The  first  data  taken  at  the  start  of  a  day's  operation  was  an  instrument 
calibration  curve.  This  was  done  by  pla'cing  calibrated  neutral  filters  of 
varying  density  serially  in  front  of  the  entrance  aperture  of  the  instrument 
and  recording  the  instrument  response  on  the  strip  chart  recorder.  After 
calibration,  the  reflectometer  was  set  to  view  the  surface  in  the  direction 
normal  to  the  surface  (0=0),  and  the  standard  smoked  magnesium  oxide 
plate  was  superimposed  a  few  millimeters  above  the  sample.  After  a  record 
of  the  instrument  output  was  obtained  for  the  standard,  the  plate  was  removed 
and  a  similar  measurement  of  the  reflectance  of  the  surface  at  0  =  0°  was 
made.  Thus  the  sample  surface  itself  as  viewed  from  6  =  0°  became  a 
secondary  standard  of  known  reflectance,  and  the  later  measurements  for 
this  particular  sample,  wavelength,  and  source  position  were  normalized 
to  the  secondary  standard  value  of  reflectance  in  the  data  reduction  process. 

After  the  above  steps  had  been  completed,  the  series  of  measurements 
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at  the  azimuth  180°  was  made  varying  the  angle  of  view  from  0  =  0°  to 
6  ^  70°  in  ten  degree  steps  and  from  0  =  70°  to  0  =  30°  in  five  degree  steps. 
When  this  was  accomplished  the  reflectometer  was  returned  to  6  =  0°  and 
the  measurement  on  the  secondary  standard  (sample  viewed  at  6  =  0°)  was 
repeated.  Thus  a  record  was  obtained  from  the  secondary  standard  once 
every  ten  minutes  or  so,  by  which  the  drift  of  the  instrument  response  or 
light  source  output  could  be  compensated  for  in  the  data  reduction  process. 

The  reflectometer  was  then  rotated  to  a  new  azimuth  and  another  scries  of 
measurements  from  6=  0°to  8  =  80°  was  begun. 

Finally,  at  the  end  of  the  series  of  sweeps  at  various  azimuths  another 
set  of  measurements  was  made  on  the  magnesium  oxide  standard  and  secondary 
standard,  after  wh' ch  either  the  wavelength,  source  position,  or  sample  was 
changed  and  the  process  was  begun  once  more  under  the  changed  conditions. 

C.  Data  obtained  from  the  measurements 

Measurements  of  the  directional  reflectance  and  degree  of  polarization 
of  the  reflected  radiation  were  made  on  each  of  ten  different  types  of  surfaces. 
The  values  of  the  various  parameters  for  which  data  were  obtained  are 
listed  as  Table  I. 

Table  I:  Summary  of  values  of  the  parameters 

for  which  reflectance  and  polarization  measurements 

wer  made  on  each  of  ten  selected  surfaces 

A  (A)  0o(°)  °> 

•4920  0,  53.  1?  78.  5  0, 180 

6430  0,53.1,  78.5  0,45,90,135,180 

7960  0,  53.  1,  78.  5  0,  180 

10250  0,  53.  1,  78.  5  0, 180 
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The  data  for  the  principal  plane  are  presented  in  the  form  of  plots  of 
directional  reflectance  versus  nadir  angle  of  observation  and  degree  of 
polarization  versus  nadir  angle  of  observation.  These  plots  make  a  total 
of  240  separate  curves  for  the  principal  plane  alone.  In  addition,  two 
hemispheric  maps,  one  for  reflectance  and  one  for  degree  of  polarization, 
are  given  for  each  of  the  ten  surfaces,  thereby  making  a  total  of  twenty 
such  hemispheric  maps.  Thus  eight  different  figures  are  given  for  each 
surface,  making  a  total  of  eighty  different  figures  for  the  measured  results. 
They  are  separated  according  to  surface  and  presented  below. 

1,  Black  loam  soil 

The  sample  of  black  loam  soil  was  obtained  in  Southwestern  Iowa  near 
the  town  of  Mt.  Ayr.  It  is  typical  of  the  soil  covering  large  areas  of  the 
Midwest,  particularly  that  of  the  corn  belt.  The  results  of  the  measurements 
are  shown  by  Figures  8  through  13  for  the  principal  plane  and  by  Figures 
14  and  15  for  the  hemisphere  at  ^  =  6430$,  Since  the  results  are  symmetrical 
with  respect  to  the  principal  plane  only  one  half  of  the  hemisphere  is  shown. 

The  other  half  is  simply  a  mirror  image  of  the  one  represented. 

2.  Desert  Sand 

The  sample  of  desert  sand  was  obtained  from  the  Mojave  Desert  about 
eleven  miles  northeast  of  Mojave.  It  is  interesting  that  where  the  sample 
was  collected  there  was  no  appreciable  change  of  darkness  of  the  material  with 
depth  below  the  surface.  This  is  not  the  case  for  a  sample  collected  previously 
from  a  spot  about  one  mile  northeast  of  Mojave.  At  that  time  it  was  noticed 
that  the  surface  appeared  to  be  bleached  out,  the  material  a  few  millimeters 
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Fig.  14  -  Hemispheric  map  of 
the  directional  reflectance  of 
black  loam  soil  (0.  =  53.  l°# 

X  =  6430$)  ° 


Fig.  15  -  Hemispheric  map  of 
the  degree  of  polarization  of 
radiation  reflected  from  black 
loam  soil  (8  =  5  3.  1°.  ^  =  6430$) 
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Nadir  Angle  of  Observation  (°) 
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Fig.  23  -  Hemispheric  map  of  the 
degree  of  polarization  of  radiation 


reflected  from  desert  sand 

(6  =  53.1°,  X  =  6430$) 


below  the  surface  being  considerably  darker  than  that  just  at  the  surface. 
The  darkness  was  not  cue  only  to  a  change  of  moisture,  as  seen  by  an 


inspection  of  the  collection  site  a  few  i  no  nth  s  later.  i  he  disturbed  area 


i  1  ■  *  _  % 


still  showed  up  as  considerably  darker  than  its  surroundings  in  that  case. 

The  desert  surface  is  obviously  different  in  the  two  locations. 

3.  Green  Grass 

The  sample  of  grass  for  which  the  data  were  taken  is  a  lush,  closely 
clipped  blue  grass  obtained  from  a  lawn  near  Philadelphia. 

4.  Dry  Grass 

This  sample  is  also  lush  bluegrass,  but  it  was  permitted  to  die  by 
lack  of  water.  There  are  a  few  blades  in  which  the  chlorophyll  has  not 
bleached  out  although  the  entire  sample  is  completely  dry  and  shows  no 
sign  of  life  at  present. 

It  can  be  seen  from  the  curves  that  there  is  a  discontinuity  in  the 
reflectance  and  polarization  at  9  =  0°.  This  comes  about  because  of  a 
preferred  orientation  of  the  blades  of  grass.  The  sample  was  taken  from 
an  operation  in  which  strips  of  sod  were  being  put  on  a  lawn,  so  it  had  been 
part  of  a  roll  of  sod.  It  can  be  seen  by  close  inspection  that  the  blades  are 
more  frequently  bent  in  the  direction  toward  the  photometer  as  the  sample 
was  installed  than  in  the  opposite  direction,  Asa  consequence  the  incident 
light  penetrated  more  deeply  into  the  grass  for  the  measurements  at  <p  =  180° 
than  at  <P  =  0°.  One  can  conceive  of  conditions  in  which  such  an  effect  would 


be  of  interest. 
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Fig.  31  -  Hemispheric  map  of 
the  degree  of  polarization  of 
radiation  reflected  from 


lipped  green  grass  ( dQ  =  53.  : 
=  6 4 30  A) 
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Fig.  38  -  Hemispheric  map  of 


the  directional  reflectance  of 
dry  grass  ( £ L  =  53.1°,  ^  =  6430^) 


Fig.  39  -  Hemispheric  map  of  the 
degree  of  polarization  of  radiation 
reflected  from  dry  grass 
<6o  =  5  3.  1°,  A  =  64  30A) 


5.  Dead  Leaves 


The  leaves  from  which  this  sample  was  prepared  were  collected  during 
the  month  of  April  from  a  wooded  site  in  VaLey  Forge  Park.  They  are 
leaves  from  the  previous  year's  growth,  and  are  mainly  from  oak  trees. 

In  order  to  make  a  reasonably  uniform  sample  as  seen  with  a  field  of  view 
covering  a  small  area,  it  was  necessary  to  break  the  individual  leaves  up 
somewhat  during  the  sample  preparation.  Otherwise  the  leaves  were  in 
their  natural  dry  state  during  the  measurements. 

6.  Broad-leafed  Plants 

This  sample  is  composed  of  small  individual  plants  of  Japanese  Spurge. 
The  leaves  are  sufficiently  dense  so  that  no  soil  is  visible  through  the 
canopy.  The  leaves  themselves  are  of  an  average  size  of  roughly  two 
square  inches,  and  do  not  appear  either  waxy  or  covered  with  fuzz.  Their 
natural  non-homogeneity  when  viewed  with  this  small  field  instrument 
causes  a  considerable  scatter  in  the  measured  data. 

7.  Crushed  Limestone  Gravel 

This  material  is  typical  of  the  foundation  often  applied  to  roadbeds. 

It  is  graded  as  3/4"  size,  although  the  individual  particles  vary  in  dimen¬ 
sions  from  roughly  1/2"  to  greater  than  3/4".  The  particles  are  generally 
angular  in  shape  as  a  result  of  the  crushing  process  by  which  they  are 
obtained  from  the  massive  limestone  rock.  This  sample  was  obtained  from 
the  quarry  located  two  miles  northeast  of  Paoli,  Pennsylvania. 
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Fig,  42  -  Directional  reflectance 

of  dead  leaves  (principal  plane, 

6  =  53.  1°) 

o 
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Fig.  47  -  Hemispheric  map  of  the 
degree  of  polarization  of  radiation 


reflected  from  dead  leaves 
(6n  =  5  3.  1°,  X  =  64 30$) 


Fig.  49  -  Degree  of  polarization 
of  radiation  reflected  from  broad 
leafed  plants  (principal  plane. 
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Fig.  52  -  Directional  reflectance 
of  broad-leafed  plants  (results  independent 
of  azimuth,  0  =  0°) 


Fig.  53  -  Degree  of  polarization 
of  radiation  reflected  from  broad 
leafed  plants  (results  independent 
of  azimuth.  0  =  0  ) 
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The  sample  of  blacktop  used  for  these  measurements  was  obtained, 
through  the  cooperation  of  the  Pennsylvania  Highway  Department,  from  a 
reconstruction  operation  on  Paper  Mill  Road  in  Oreland,  a  suburb  of 
Philadelphia.  The  road  has  fairly  heavy  traffic  but  since  there  was  no 
painted  centerline  there  were  no  preferred  strips  of  tire  wear  along  the 
roadway.  The  sample  surface  is  not  particularly  smooth,  asphalt  covered 
pebbles  protruding  as  much  as  one  eighth  to  three  sixteenths  inches  above 
the  general  level  of  the  surface.  There  was  no  definite  information  avail¬ 
able  on  the  length  of  time  the  road  surface  had  been  in  place,  but  it  was 
"several  years" . 

9  Concrete 

The  sample  of  concrete  for  which  the  measurements  were  made  was 
obtained  from  Wayne,  Pennsylvania.  It  is  a  section  of  well-weathered 
sidewalk.  The  general  surface  is  relatively  smooth  as  concrete  surfaces 
run,  and  some  small  pebbles  of  the  order  of  one  eighth  inch  dimater  are 
visible  in  the  matrix. 

10.  Canvas 

The  sample  of  canvas  used  is  that  supplied  by  the  Naval  Ordnance 
Test  Station,  China  Lake,  California.  It  appears  to  be  a  weathered 
piece  of  a  standard  type  of  green-colored  canvas. 
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Fig.  70  -  Hemispheric  map  of 
the  directional  reflectance  of 
weathered  blacktop  (0Q  =  53,1°, 
X  =  6430&) 


Fig.  71  -  Hemispheric  map  of 
the  degree  of  polarization  of 
radiation  reflected  from  weathered 
blacktop  (6  =  53.1°,  X  =  6430A) 
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Fig.  79  -  Hemispheric  map  of  the 
degree  of  polarization  of  radiation 
reflected  from  weathered  concrete, 
(6  =  53.  i°,  X  =  6430A) 
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Fig.  80  -  Directional  reflectance 
of  canvas  (principal  plane, 

6  «  78.  5°) 
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IV.  Theory 

The  aim  of  the  following  sections  is  to  incorporate  the  laboratory 
measurements  of  the  directional  reflectance  and  degree  of  polarization 
of  the  reflected  radiation  for  a  given  sample  into  the  solution  of  the 
radiative  transfer  equation.  This  enables  one,  with  reasonable  approxima¬ 
tion  to  compute  the  intensity  and  the  degree  of  polarization  of  the  radiation 
emerging  from  the  top  of  the  atmosphere.  Finally,  the  effects  of  the 
atmosphere  on  the  contrast  inherent  between  an  object  and  the  surrounding 
background  will  be  shown.  Before  proceeding,  it  seems  wise  to  mention 
the  assumptions  underlying  some  of  the  above  statements,  and  to  review 
the  work  on  contrasty  due  to  R.  Fraser  (1964). 

In  the  late  1940's,  S.  Chandrasekhar  was  able  to  formally  solve  the 
equation  of  radiative  transfer,  including  the  effects  of  polarization  and 
multiple  scattering.  The  solution  appeared  as  a  series  of  articles  in  the 
A  strophysical  Journal  which  were  collectively  presented  in  his  treatise  of 
1950  titled  Radiative  Transfer.  The  most  important  restrictive  assumption 
to  realize  is  that  the  solution  obtained  has  relevance  for  a  so-called  Rayleigh 
atmosphere,  i.  e,  one  in  which  the  scattering  particles  are  very  small 
compared  to  the  wavelength  of  radiation.  However,  a  "real"  atmosphere 
contains  larger  particles,  usually  referred  to  as  aerosols,  and  composed  of 
duet,  by  -products  of  combustion  processes,  etc.  .  Attempts  to  include  the 
presence  of  aerosols  in  the  equation  of  radiative  transfer  have  teen  made  by 
several  authors,  in  particular  Z.  Sekera  (1956),  R.  Fraser  (1959),  and 
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D.  Deirmendjian  {195?;.  The  problem  is  complicated  by  several  practical 
considerations.  In  the  first  place,  one  must  have  good  measurements  of 
the  size  distribution  of  aerosol  particles  per  rnit  volume.  This  imm¬ 
ediately  raises  the  question  of  whether  this  distribution  remains  invariant 
with  altitude,  as  well  as  the  nature  of  the  vertical  density  profile  of 
aerosol  content.  It  was  also  found  that  the  aeiosol  content  of  the 
atmosphere  is  dependent  on  geographical  location,  and  even  at  the  same 
locality  may  vary  with  time.  The  approximations  usually  resorted  to, 
which  will  be  followed  in  this  report,  assume  that  the  density  of  aerosol 
particles  decreases  exponentially  with  altitude,  and  that  the  size  distribu¬ 
tion  of  aerosol  particles  per  unit  volume  is  constant.  These  assumptions 
simplify  the  transfer  equation  considerably  once  ar.  aerosol  model  is 
chosen.  We  will  elaborate  on  this  later  in  the  text. 

Another  important  assumption  introduced  by  Chandrasekhar  in  order 
to  make  the  problem  more  tractable  is  that  the  atmosphere  is  idealized  in 
that  it  is  considered  to  be  plane  parallel.  In  other  words,  the  atmosphere 
is  considered  to  be  bounded  by  two  parallel  planes  of  infinite  extent  such 
that  the  characteristics  of  the  atmosphere  depend  only  on  the  altitude. 

This  idealization  avoids  the  problem  of  sphericity  of  the  atmosphere  and 
Earth,  and  the  assumption  is  valid  as  long  as  one  is  not  interested  in 
atmospheric  effects  near  the  horizon. 

As  an  application  of  his  solution  to  the  molecular  scattering  problem, 
which  was  originally  solved  assuming  that  there  was  no  ground  present, 
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(that  is,  the  atmosphere  extended  indefinitely  downward),  Chandrasekhar 


introduced  the  effects  that  the  presence  of  a  ground  would  have  on  the 
solution  that  he  obtained.  To  be  more  specific,  a  Lambert  surface  was 
chosen  for  two  reasons.  First,  since  a  Lambert  surface  is  one  which 
reflects  incident  radiation  with  the  same  intensity  in  all  directions,  and 
in  addition,  the  reflected  light  is  unpolarized,  the  resulting  equations 
were  very  much  simplified.  Secondly,  it  was  believed  that  the  approxima¬ 
tion  of  a  real  surface  by  a  Lambert  surface  was  sufficiently  accurate  for 
most  practical  purposes.  It  was  shown  that  the  presence  of  the  ground 
created  an  additive  correction  to  the  so-called  "standard"  solution; 
that  is  the  solution  obtained  without  the  presence  of  a  ground.  In  I960, 
extensive  tables  were  published  (Coulson,  Dave,  and  Sekera)  giving  the 
numerical  solution  of  the  transfer  problem  with  a  Lambert  surface 
according  to  the  theory  oi’  Chandrasekhar.  Fraser  has  uced  these,  tables 
in  determining  the  atmospheric  effects  on  contrasts.  Our  major 
innovation,  of  course,  will  be  to  introduce  a  non- Lambertian  surface. 

Before  doing  this,  it  seems  appropriate  at  this  time  to  present  a 
statement  of  the  problem,  and  to  .roduce  the  necessary  notation.. 

In  figure  88,  we  have  a  simplified  representation  of  the  physical 
aspects  of  the  problem.  A  small  target,  resting  on  an  infinis*.,  homogeneous 
background,  both  illuminated  by  sunlight,  is  viewed  by  an  observer  (cr 
instrument)  situated  at  the  top  of  the  atmosphere.  The  field  of  view  is 
larger  than  the  solid  angle  subtended  by  the  target,  and  therefore  a 


109 


OBSERVER 


,V,?  //  Q  Q  a  H£  o  vS  fOV'VO 

Figure  88  «  A  Simplified.  Representation 
**C  l Re  Cowtradt  Determination  Problem 


"contrast”  between  the  target  and  the  background  is  perceived.  It  is 
natural  to  call  this  an  apparent  contrast,  and  reserve  the  term  intrinsic 
contrast  for  that  contrast  an  observer  would  perceive  at  the  bottom  of  the 
atmosphere.  In  other  words,  the  intrinsic  contrast  is  altered  due  to  the 
presence  of  the  atmosphere. 

The  notation  to  be  employed  in  what  follows  is  listed  below.  Detailed 
explanation  of  some  of  these  items  will  be  postponed  until  the  need  for 
their  introduction  arises. 


0O  =  the  zenith  angle  of  the  sun  (as  seen  from  the  target) 
or  the  nadir  angle  of  the  target  as  seen  from  the 
direction  of  the  sun 

2  cos  eo 

6  s  the  zenith  angle  of  the  observer  (as  seen  from  the  target) 

or  the  nadir  angle  of  the  target  (as  seen  from  the  observer) 

p  =  cos  6 


<fi  =  the  azimuth  or  horizontal  angula  r  measi  re 

of  the  observer  relative  to  the  sun  (i.  e.  the 
sun  is  at  azimuth  <p  =  0°).  In  the  figure,  the 
observer  is  at  <p  -  180°  locking  in  the  direction 
<p  =  0°. 

z  =  altitude 


X  = 

8  (X,z) 


A 

B 


wavelength  of  radiation 

the  volume  scattering  coefficient  for  the  Rayleigh 
or  "molecular"  part  of  the  atmosphere.  This  is 
the  rate  at  which  incident  energy  is  scattered  per 
unit  volume  per  unit  path  length.  It  depends  on 
wavelength  and  density,  which  in  turn  implies 
an  altitude  dependence. 

the  corresponding  volume  scattering  coefficient 
for  the  aerosol  content  of  the  atmosphere. 


R 

T  =  the  normal  optical  thickness  of  the  molecular 

component  of  the  atmosphere. 

R  °°r  R 

By  definition,  T  =  J  (3  (z)  d  z, 

o 

A 

T  =  the  normal  optical  thickness  of  the  aerosol 

component  of  the  atmosphere. 

A  ”  A 

T  =  J  /9  (z)  d  z 

o 

Tj  s  the  total  normal  optical  thickness  of  the 

atmosphere. 

R  A  “  R  r  A 

T  *  T  +  T  =  J  /S  (z)dz+J  /3  (z)  d  z 

o  o 


Ij.  (0;  »  (P)  5  tne  emergent  intensity  of  the  target-reflected 

radiation  at  the  top  of  the  atmosphere  in  the 
direction  p  .  <£• 


<p)  =  the  intensity  of  the  target  reflected  radiation 

at  the  ground 


Ifi{0 


the  corresponding  intensities  for  the  homogeneous 
background. 


The  quantities  p  (p  ,  <P)»  P  (p  ,  <p)  , 


have  been  introduced  previously, 
and  their  complete  description  will 
not  be  repeated  here.  However,  it 


is  important  to  remember  that  Ig  and  Ir  are  the  intensity  components 


parallel  and  perpendicular  to  the  vertical  plane  passing  through  the 


zenith  and  the  azimuth  of  the  observer.  Later  on,  we  will  have  need 


for  intensity  components  parallel  and  perpendicular  to  the  plane  of 


polarization*  The  notation  T  and  lj  will  be  reserved  for  this. 

y  s  the  scattering  angle.  This  is  the  angle  formed 

by  the  incident  and  reflected  beams.  The  plane 
irmined  by  these  beams  is  referred  to  as  the 
scattering  plane. 

_A 

L  (7)  These  are  normalized  matrix  elements 

-  describing  the  scattering  properties  of 

iJ*  (y)  a  particular  aerosol  model. 

1 2 

C  (0;p,<p)5  apparent  contrast  at  the  top  of  the  atmosphere, 
and  in  the  direction  Ji, 

C  (Xj ;  p,  <p)  £  intrinsic  contrast  at  the  bottom  in  the 
direction  p  ,  <p  . 

y  (Tpp,  <p)  =  the  contrast  attenuation  coefficient,  which 
holds  in  the  absence  of  an  analyzer  in  the 
observer's  optical  system. 

ye(Ti,  ,  <p)=  the  contrast  attenuation  coefficient  which 
holds  if  the  transmission  plane  of  an 
analyzer  in  the  observer's  system  is 
oriented  perpendicular  to  the  plane  of 
polarization. 

yr(Ti;M»  5  t*ae  corresponding  contrast  attenuation 

coefficient  which  holds  if  the  transmission 
plane  of  an  analyzer  is  oriented  parallel  to 
the  plane  of  polarization. 

Having  introduced  the  notations  for  contrast  and  the  contrast 
attenuation  coefficient,  it  remains  to  define  these  terms  and  point  out 
the  great  utility  of  the  concept  of  the  contrast  attenuation  coefficient.  The 
following  discussion  of  these  parameters  is  patterned  after  that  of 
Fraser  (1964).  The  contrast  at  any  altitude  z  corresponding  to  a  normal 


optical  depth  T,  and  in  the  direction  is  defined  as: 


1 


C(T;M,<P)  5  I t(r-,n,<p)  -  I b(t;m.«P) 

i" 

In  particular,  the  intrinsic  contrast  is: 

G(Ti;M.<P)  =  -  IB  (lyM.P) 

Ifi  (Tj  ;  M.  <P) 

The  apparent  contrast  as  perceived  at  the  top  of  the  atmosphere  is: 

C  (0 -,H,<p)=  It  (0 ;M,<p)  -  IB(0;  M.<P) 

1  (0;  H,  <p) 

B 

The  emergent  intensity  I  (0;  pi,  (p)  is  related  to  the  intensity  reflected 
from  the  ground,  I  (r^  pi,  <p),  by  integrating  the  transfer  equation.  The 
result  may  formally  be  written: 

-Tl  /M  T1  -(/M 

(1)  I  (0;  p)  =  I  < Tj ; M  .  <f> )  e  +  J  J(t;  P* . <P>  e  dt 

o  M 

The  quantity  J  is  usually  referred  to  as  the  source  function,  and  the 
integral  term  is  responsible  for  multiple  scattering  effects.  The  as¬ 
sumption  of  very  small  target  size  implies  that  the  contribution  of  the 
target  to  multiple  scattering  effects  is  minimal.  Hence  if  one  assumes 
that  the  multiple  scattering  effects  above  the  target  are  the  same  as  above 
the  surrounding  background,  then  the  integral  terms  for  both  the  target 
and  background  will  be  the  same.  Therefore,  substituting  equation  (1) 
into  the  definition  of  apparent  contrast  at  the  top  of  the  atmosphere,  we 


have: 


(2)  C(0;M,<P)  =  |  It  (TpM.V)  -  Ib(ti;M.V)} 


-  Tl/, 


IB(ri;^,<p)  e  lyV  +  j*  j  (t;#i.V>)  e  M  dt 

o  M 

The  contrast  attenuation  coefficient,  y  (Ty,p,<p)  is  now  defined  such 

that: 

(3)  C(0;m.<P)  =  y  (tl;M,P)  C  (Ti;/l,v) 

In  other  words,  the  apparent  contrast  is  the  product  of  the  intrinsic 

contrast  and  the  contrast  attenuation  coefficient.  Hence  y  (7 \  ;  fx  ,ifi)  is  a 

measure  of  the  reduction  of  contrast  due  to  atmospheric  effects.  Since 

C  (Ti;#*lV)  =  I  (Ti;m,<P)  -  Ir^M.P)  .  ,  .  _ 

*  1  _ _ 2 _ ,  have  by  substitution, 

ib(ti;  <p) 

■  ri/ 

y(ri;M,<p)  =  C  (0;m.Q)  =  e  M 

C  d>) - - T - ~~t/ 

IB  <T1  e  1/m  +  /  J  (t;M»  <P)  e  Mch 


The  utility  of  y  can  be  seen  from  the  fact  that  the  resulting 

expression  is  independent  of  the  target,  and  depends  only  on  the  back¬ 
ground.  The  characteristics  of  the  target  are,  of  course,  introduced  in 
the  calculation  of  the  intrinsic  contrast. 

It  is  further  seen  from  the  expression  for  the  contrast  attenuation 
coefficient  that  0  £y  (Tp  s  1.  If  the  atmosphere  is  optically  thick 

'Vu 

so  that  Tj  is  relatively  large,  then  (T^;  gf,V>)  e  'will  be  small, 

Ii  -t/M 

while  J  (t;  p,C)  e  dt  may  become  quite  large,  giving  rise  to  a 

o  ~ji~ 

small  contrast  attenuation  coefficient.  In  this  case,  the  intrinsic  contrast 
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is  greatly  reduced  as  th<_  light  traverses  it  s  path  through  the  atmosphere. 
On  the  other  hand,  if  the  atmosphere  is  optically  thin,  so  that  T.  is 


relatively  small,  then  1,.  (t,  c  will  be  large  compared  to 

V  -t  /,  b  1 

J  J  (t;  p.O)  e  dt  ,  at. a  in  this  case,  we  have  a  large  contrast 


attenuation  coefficient.  Hence  the  intrinsic  contrast  is  not  greatly 
reduced  as  the  light  passes  out  through  the  atmosphere.  Of  course,  these 
qualitative  descriptions  are  what  one  would  intuitively  expect  in  an 
optically  thick  or  optically  thin  atmosphere. 

At  this  point,  it  is  useful  to  distinguish  between  intensity  of  radiation 
which  is  directly  transmitted  and  that  which  is  diffusely  transmitted.  A 
convenient  notation  for  this  purpose  is  described  as  follows:  In  general, 
the  emergent  intensity  may  be  considered  to  be  composed  of  five  separate 
components.  The  word  "component"  here  is  used  synonymously  wi*h 
"contribution",  and  should  not  be  confused  with  the  previous  reference  to 
vector  components  of  intensity. 

I  =  the  solution  to  the  radiative  transfer  problem  in  the  absence  of  a 
"ground".  This  is  the  so-called  standard  solution  alluded  to 
previously . 

I  55  that  contribution  to  the  emergent  radiation  which  consists  of  light 
which  was  transmitted  directly  downward  through  the  atmosphere, 
ref’ected,  and  transmitted  directly  outward  through  the  atmosphere. 

I  “  that  contribution  to  the  emergent  radiation  which  consists  of 
light  which  was  transmitted  diffusely  downward  through  the 
atmosphere,  and  transmitted  directly  outward  through  the  atmos¬ 
phere. 

I  =  that  contribution  to  the  emergent  radiation  which  consists  of  light 
which  was  transmitted  directly  downward  through  the  atmosphere 
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and  transmitted  diftuscly  outward  through  the  atmosphere. 

^dd  ~  that  contribution  to  the  emergent  radiation  which  consists  of 

light  which  w-as  transmitted  diffusely  downward  through  the 
atmosphere,  and  transmitted  diffusely  outward  through  the 
atmosphere. 

Before  discussing  these  components  in  detail,  and  presenting  the 
expressions  which  will  be  used  for  each,  it  is  worthwhile  to  note  that  the 
contrast  attenuation  coefficient  can  now  be  written  as: 

y(Ti;|i,<p)  =  4  :dD 

V".*)  +  W*‘-<?)+  1ddM,(S)  +  lD^V)  1 

(For  ye  (T^  ;(i,<p)  and  yr(Tj  we  need  compute  the  i  and  j  components 

of  each  of  the  above  five  "contributions".  } 

As  a  first  step  in  computing  these  intensitities,  it  is  necessary  to 
determine  the  normal  optical  thickness  which  applies.  For  the  molecular 
component  of  the  atmosphere,  the  following  table  is  taken  from 
Deirmendjian  (1955): 

Table  II:  The  normal  optical  thickness,  T, 

as  a  function  of  wavelength,  X  . 

T  1.00  0.50  0,25  0.15  0.10  0.05  0.02 

X  (A)  3120  3715  4365  4950  5460  6440  8090 

In  order  to  illustrate  the  computation  of  y  (Tt  ;ft,<p),  we  have  chosen 

a  particular  sample  More  specifically,  we  shall  introduce  the  surface 

reflection  and  polarization  characteristics  of  desert  3oil  for  the  wavelength 

o  o 

X  =  6430A,  and  a  zenith  sun  angle  of  6^  =  53.1  (MQ  =  .60).  Hence, 
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for  this  wavelength,  the  required  molecular  normal  optical  depth  is 


A 

In  determining  T  ,  the  normal  optical  depth  of  the  aerosol  content, 
one  must  choose  a  model  for  the  size  distribution  of  aerosols  per  unit 
volume.  We  have  chosen  a  ''continental"  model  aerosol  from  .'ime-rig 
several  considered  by  Fraser  (1959)  as  being  most  representative  for  our 
purposes.  It  is  characterized  by  the  following  size  distribution: 

Interval  of  radius,  a  Nun^ber  of  particles  Number  of 

(micron)  per  cm  of  air  per  micron  particles  per  cm^ 


of  radius 

Lower 

Upper 

4 

.  03 

0.  1 

2  251  x  10 

-4 

1575.  7 

0.  1 

20.  0 

2.  251  x  a 

753.  3 
2329-  0 

As  mentioned  previously,  this  surface  value  will  be  assumed  to  decrease 
exponentially.  For  this  model,  Fraser  computed  the  volume  scattering 
coefficients  at  the  surface  for  three  different  wavelengths.  His  result 

.  c  A  -6-1 

for  =  6250A  isP  =  1.  06  x  10  cm  .  By  interpolating  according  to  an 

inverse  wavelength  dependence,  we  have  for  ^  =  6430A,  that  P I (0)  = 

-6  -1 

1.03  x  10  cm  .  Finally,  utilizing  the  assumption  of  an  exponential 

A  A  -z/H 

density  decrease  with  altitude,  we  have  that  P  (z)  =  /?(0)  e  .  The 
parameter  H  is  called  the  scale  height.  Physically,  it  is  that  alticude  at 
which  the  volume  scattering  coefficient  assumes  l/e  of  its  surface  value. 
The  value  H  =  .98  kilometers  is  given  by  Penndorf  (19?4). 
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A 

By  definition,  T 


e 


"  A  “  A  -Z/H  A 

-  J  /3  (z)  dz  =  J  /3  (0)  e  =  0  (0;  H  r. 

O  o 


-6  -i  3  .  ,  Q 

1,  03  x  10  cm  x  .98  x  10  cm  =  .101  .  Therefore,  for  .X  =  6430  A, 


-1 


the  total  normal  optical  thickucss, 


R  fi 

=  T  4  C  is  approximately  equal 


to  .  15. 


We  shall  now  discuss  briefly  each  of  the  five  components,  beginning 
with  I  (0;  Ji,  <p).  Our  aim  is  to  determine  the  I  and  I  components  for 
each  of  the  five.  To  reiterate,  the  subscript  s  is  for  "standard'',  the 
reference  being  to  the  solution  of  the  radiative  transfer  equation  in  the 
absence  of  a  "ground".  In  the  present  problem,  Ig  consists  of  contribu¬ 
tions  from  both  molecular  scattering  and  aerosol  scattering.  We  have 

made  the  simplifying  assumption  that  these  scattering  processes  act 

R  ft 

independently,  so  that  Ig  can  be  written  as  Is  =  I8  +  Is  ,  separating 

the  molecular  and  aerosol  contributions.  Hence  we  desire  the  solution 
R 

for  Ig  with  a  normal  optical  thickness  of  T  =  .  05,  and  the  solution  for 

A  K 

I  with  a  normal  optical  thickness  of  T  =  .10.  Values  of  Ig  (0;  /!,<&)  for 

T  =  .  05  are  given  in  the  tables  of  Coulson,  Dave,  and  Sekera  (I960  ). 

Actually  the  Stokes  parameters  I  and  Q,  are  given  which  make  a 

R 

determination  of  I„  and  1^  for  Is  trivial.  It  is  to  be  noted  that  this 

v  *  “ 

solution  includes  the  effects  of  multiple  scattering.  In  compiling  these 
R 

tables  for  Ig  ,  an  assumption  was  made  that  the  flux  of  parallel  radiation 
incident  on  the  top  of  the  atmosphere  is  set  at  15  Fc  (F0  =  1)  units,  per  unit 
area  oriented  normal  to  the  direction  of  propagation,  which,  of  course. 
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makes  an  angle  of  Qq  with  the  zenith.  This  assumption  will  be  followed  in 
the  present  report. 

fl 

In  computing  I  ,  we  follow  the  usual  approximation  of  neglecting  the 
effects  of  multiple  scattering,  since  their  inclusion  renders  the  problem 
intr-ctabie.  The  transfer  equation,  including  primary  scattering  may  be 
written: 


H  dl_ 

dr 


“T/u  p 
I  -  1/4  e  P 


where: 


I. 

+  I. 

= 

1 

J 

I. 

-  I  1 

„  1 

J  J 

A-  Ka>  4 ^ 


hzM  hiM 


F  =  1 

0 

A 

The  scattering  matrix  P  describes  the  aerosol  scattering  of  incident 
radiation,  and  enables  one  to  compute  the  intensity  components  Ij  and  Ij 
parallel  and  perpendicular  to  the  scattering  plane,  which  is  normal  to  the 

A  A 

plane  of  polarization.  The  elements  (y  )  and  L-j ^ ( y)  are  the  normalized 

matrix  elements  introduced  previously.  Fraser  (1959)  has  computed 

these  elements  for  the  "continental"  aerosol  model,  and  for  a  wavelength 
o 

X  =  6250  A.  Since  they  do  not  change  rapidly  with  wavelength,  these 


values  are  also  appropriate  for  X  =  6430A. 


120 


An  important  fact,  pointed  out  by  Deirmendjian  (1959),  is  that  the  assumption 

A 

of  exponential  density  decrease  of  aerosol  content  ensures  that  BL-(y)  =  0  , 

B  r 

This  means  that  in  solving  the  transfer  equation  above,  these  elements  can 
be  brought  out  from  under  an  integral  with  respect  to  T.  In  fact  since. 


A  - 

r  a  *  I 

| 

A 

P  F  = 

Lu(y)  l12  (y) 
if12(y )  Lu(y) 

i 

0 

■  ’ 

Ln  <« 
4  M 

The  transfer  equation  assumes  the  form: 

-  r/ 


- 

—  “ 

M  d 

I 

I 

_ ; 

_  Q  . 

= 

d  r 

1/4  e 


'Mo 


r  t 

l  (y) 

li 

4<y) 


or 


TJ 

>  1 1 


(5a)  fu  dl_  =  I  -  1/4  Lu  (y)  e  '"o 

dr 

-  '  T/,m„ 

(5b)  M  dQ_  =  Q  “  1/4  L  (y)  e 

dr 

The  solutions  to  these  linear  differential  equations  may  be  written, 
remembering  that  the  boundary  conditions  are  such  that  there  is  no  radiation 

4 


reflected  from  a  "ground",  as: 


(6a)  I  =  lj  +  Ij  =_1 _  Ln  (y)  £  Mo  J  1 

4M  M  +  Mo  I 


’T  {-  + 
e  M  Mo 


(6b)  Q  =  I. 


I  =  1 
J 


4M 


l  (y)  MMC 


M  +  M^ 
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I. 


I 

j 


=  1  P^O 

8  p 

=  1 

ITp  M  4  Pp 


j  Lnm  4  Lii(>) 


L  (>)- 
11 


Lu<y)1| 

t 


These  components  arc  parallel  and  pe  rpendicula  r  to  the  scattering 

plane.  To  transform  these  components  to  components  parallel  and 

perpendicular  to  the  vertical  plane  through  the  observer's  azimuth,  we 

2 


need  subject  them  to  the  matrix  transformation 
Hence  we  have. 


2 

sin-,  X 

cos  x 


cos^  X 
sin  X 


r  r  i 

.  ?  2 

T  2  2 

e  1 

sin  X  cos  X 

h 

- 

I.  sin  X  +  IjCOs  x 

i " 

2.  .  2 

,  2  ,  .  2V 

i 

i 

cos  X  sin  X 

i  i- 

I  cos  X  +  I  .sin  X 

i  r 

i 

j  . 

i 

1  .5 

L  —I 

In  the  computation  of  Ipjpj  <P  )  we  must  introduce  our  laboratory 

measurements,  P  )  and  P(  ).  If  ^Fn  units  of  flux  per  unit 

area  normal  to  the  direction  of  propagation  are  incident  on  the  top  of  the 

atmosphere  at  an  angle  relative  to  the  zenith,  then  the  directly  transmitted 

~T\/ 

radiant  flux  reaching  the  ground  per  unit  horizontal  area  is  ^  F0  e  1 '  P0  * 

The  intensity  of  the  radiation  reflected  into  a  unit  solid  angle  in  the  pi  v 

ref  “Tl/u 

direction  would  be  I  (M,<P)  =  P  (P,<P)  tt  Fq  M0  e  '  ^o  .  We  now  make 
use  of  the  definitions  of  the  Stokes  parameters: 


I  =1  +  I 

e  r 


Q  -  I 


e  r 
U  =  (I  -  Ir)  tan  2  X  , 


and  P  =  (Q2  4  U2  )1/2 


It  will  be  remembered  that  X  the  angle  between  the  plane  of  polarisa¬ 
tion  and  the  vertical  plane  through  the  observer's  azimuth.  The  assumption 
that  the  plane  of  polarization  is  normal  to  the  scattering  plane  enables  one 
to  compute  X  by  spherical  trigonometry.  Hence  the  values  of  and 

determine  the  angles  X  (and  y )  directly.  We  then  have  the  following  relation¬ 


ships: 


(  2  2  V 

P  (V,<P)  =  ]Q  (H,V)  +  U  (M,<P)) 

I 


)  2  Z 

|  Ie  (M.V)  -  lr  (>i.V>))  (le  -  Ir  tan2  2  X 


I  OV?) 


f  Ie  (M,<P)  -  Ir  (M.<P)j  sec  2  X 


I  (M,<P) 

and  P  (M,P)  I  cos  2  X  =  1-1  =  (I  -  ir)  -  I 

e  r 


=  1-21 


or  2Ir  =  I 


1  -  P  cos  2  X 


Therefore,  it  follows  that 


!(eref>  =  i/2  !ref  (^,<p)  |  1  +  p  (M,V)cos2x] 


I(/ef)  (P,«P)=l/2I(ref) 


-  P  m.«P)  cos  2x]. 

-T, 


"  1  1  / 

Each  of  these  components  is  attenuated  by  the  factor  e  P  in  its  out¬ 
ward  passage  through  the  atmosphere.  Substituting  the  expression  for 
ref 

1  (M.VJ.we  have  the  final  expressions: 
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(7a)  ^Ie(P.V>)j  =  1/2  P  (p  ,P)  TTp0M0  e  l/P0  Vp  |  1  +  P  (P.V5)  cos  2  \J 

_7  _T .  | 

(7b)  lr  (P.pJdd  =  1/2  P(M.P)»F0P0e  1/M°  e  ^  j  l  -  P  (P,<P)  cos  2  \J  . 

In  considering  the  remaining  three  components,  we  have  made  the 
assumption  that  the  scattering  processes  for  these  components  obey  the 
Rayleigh  or  molecular  scattering  law.  This  approximation  is  probably  not 
a  bad  one  for  two  reasons.  First,  the  importance  of  these  components 
diminishes  with  decreasing  optical  thickness,  and  secondly,  the  aerosol 
scattering  matrix  approa-^es  the  Rayleigh  scattering  matrix  as  the  particle 
size  approaches  zero.  The  effect  of  the  aerosol  content  is  to  increase  the 
optical  depth  to  =  .  15.  The  expres  sion  for  I  (0 ;  p  ,  <p)  is 

— »  V  rr  '  '  (ref)  *  '  '  ' 

(8)  I  (0  ;P,P)  =  1  J  J  I  (T;p,<p;.U  ,<p)  I  (P.tP)  dp  d<P 

Dd  W  o  o  ~ 

There  is  little  need  in  describing  the  computational  schemes  involved 

in  this  integral,  as  the  details  are  available  in  an  appendix  to  General 

Electric  Technical  Information  Series  R64SD74  (1964).  The  important 

( l*  *f ) 

thing  to  realize  is  that  the  term  F  e  contains  the  measured  values  of  P 
and  P,  which  are  introduced  into  the  intensity  components  in  the  manner 

described  in  the  section  on  (P .<£). 

-  -  uu  “ 

For  the  remaining  two  components,  I  and  I  ,  one  further  assumption 

dD  do 

was  imposed.  Because  of  the  difficulty  in  evaluating  the  downward  diffusely 

transmitted  light  in  the  presence  of  a  non-Lambertian  surface,  it  was 

decided  to  idealize  the  surface  to  a  Lambert  surface  reflecting  with  an 

1  rr 

albedo  of  IT  =  J  J  P  (P,  <P )  P  d  p  d  <P  .  This  approximation  is  a  good  one 
o  o 
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for  two  reasons.  First,  as  noted  before,  these  components  are  in  general, 

small.  Secondly,  a  computation  of  the  downward  diffusely  transmitted 

light  would  involve  an  integration  over  the  whole  hemisphere,  thus 

originating  a  tendency  toward  averaging.  Actually,  the  problem  is  much 

more  difficult  than  it  appears,  since  an  exact  solution  would  require  the 

solution  of  certain  integral  equations.  In  summary  then,  the  components 

I  and  1  will  be  approximated  by  introducing  a  Lambert  surface  with 
dD  dd 


albedo  R.  These  intensities  can  readily  be  written  in  tc  .ms  of  the 
function  >'e  ,  y  ,  and  "5  introduced  by  Chandrasekhar  in  the  solution  of 
the  Lambert  surface.  For  the  sake  of  completeness,  these  expressions 
are  given  as  follows: 


(9)  I  <«  =  u  v  Pc  )  y  <m0i  +  y r<M0> 


2<!'r  s' 


~T‘/M  1  ^  *T,/» 

e  o  r  R  e 

J  7T 


<'°i  >dd  (m)  =  J5.M  "F0  f  »  >,  (n.)  _  C"T'V0 


Te(^)  +  ?r(M)  -  2e 


•r,/. 


The  functions  and  s  have  been  tabulated  by  Sekera  and  col- 

e  r 

laborators  (1952  )  for  T  =  .15.  Two  things  are  noteworthy  about  the  above 

expressions.  First,  they  are  azimuth  independent.  Secondly,  their 

polarization  is  almost  negligible  so  that  the  I  and  I  components  of  each 

e  r 

can  both  be  taken  equal  to  one  half  the  value  of  the  intensity. 

This  concludes  a  brief  description  of  the  various  components  considered 
in  the  calculation  of  the  emergent  intensity.  Asa  matter  of  interest,  the 
relative  contribution  of  these  various  components  in  the  sun's  vertical  is 
presented  as  figure  89-  With  these  data,  one.  can  now  compute  the  contrast 
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attenuation  coefficient,  as  previously  described. 

One  further  consideration  remains,  namely  the  computation  of  yfi 
and  y^  at  other  azimuths.  The  intensity  components  Ie  and  Ir  have  been 
computed  relative  to  the  vertical  plane  through  the  observer's  azimuth. 
What  we  now  require  are  the  intensity  components  parallel  and  normal  to 
the  plane  of  polarization.  If  the  observer  is  situated  in  the  sun's  vertical 
plane,  as  is  the  case  in  the  preceding  figure,  then  there  is  no  complication, 
since  the  sun's  vertical  coincides  with  the  scattering  plane  and  is  normal 
to  the  plane  cf  polarization.  If  this  is  not  the  case,  then,  in  terms  of 
angle  X,  we  have 

( s)  2  2 

(11a)  I.  =  I  =  I  cos  X  -  Ir  sin  X 

j  r  ® _ * _ 

cos^  X  -  sin2  X 


(s)  2  2 

(lib)  I.  =  \  =  Ircos  X-  Iesin  X 

cos^X  -  sin^X 


iiiH< 


vf  r*  Arvtnnf  ' 


1 1  a  r»  « 


The  results  of  the  computations  for  the  contrast  attenuation  coefficients 

are  displayed  in  figures  9Q-94.  As  mentioned  previously,  these  computations 

are  for  the  "Desert  Soil"  sample,  for  a  wavelength  ^  =  6430^,  and  for  a 

o 

sun  zenith  angle  M0=  .  60  (Bq  =  5  3. 1),  Figure  90  shows  the  values  of  ye, 
y,  and  yr  in  the  principal  plane  (i.  e.  the  vertical  plane  containing  the  sun 
and  the  zenith).  A  point  of  interest  is  the  fact  that  maximum  values  for 
both  and  y  are  achieved  on  the  solar  side  of  the  nadir.  Figures  91  ,  92, 
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Figure  91  -  The  contrast  attenuation 
coefficient  yr  as  a  function  of  azimuth, 
for  Desert  Soil,  0  =  53fl,  and  ^  =  6430A 


NADIR  ANGLE  (°) 


Figure  93  -  The  contrast  attenuation 
coefficient  y  as  a  function  of  azimuth, 
for  Desert  So  l,  6  =  53?1  and  6430  % 


and  93  are  an  extension  of  figure  90  and  show  the  azimuth  dependence  of 
V  v  and  v_  re snertivelv.  Aoain.  the  v  and  V  plots  clearly  show  the 

-  d  '  - ^  g  i  J  O  ~m  1  <  g  *  * 

asymmetry  with  respect  to  the  90°  azimuth. 

While  these  figures  show  the  direction  of  maximum  contrast  attenua¬ 
tion  coefficient,  it  cannot  be  inferred  that  the  apparent  contrast  is  maximum 
in  this  direction,  since  the  apparent  contrast  is  the  product  of  the  "intrinsic" 
contrast  and  the  contrast  attenuation  coefficient.  However,  in  the  absence 
of  any  knowledge  concerning  the  "intrinsic"  contrast,  the  curves  do  show  in 
■what  directions  the  reduction  in  contrast  due  to  atmospheric  effects  can  be 
minimized. 

Finally,  the  relative  advantage  of  using  a  polarizer  is  given  as  a 

function  of  direction  in  figure  94.  The  relative  advantage  is  defined  by 

ye  -  y  x  100.  A  very  strong  gradient  is  seen  toward  the 

y  (M  ,<P) 

horizon  in  the  90°  azimuth  direction.  However,  as  noted  in  the  previous 
figures,  the  contrast  attenuation  coefficients  approach  zero  as  one  looks 
toward  the  horizon. 

VI.  Recommended  Further  Research 

As  can  be  seen  from  the  above  results  and  discussion,  a  good  start 
on  an  investigation  of  the  problems  of  object  detection  by  polarized  light 
has  been  made  on  this  project.  There  are,  however,  many  aspects  to  the 
task,  and  it  has  not  been  possible  within  the  scope  of  the  present  effort  to 
cover  all  of  the  areas  adequately.  In  view  of  the  importance  of  this  new 
and  relatively  unexploited  approach  to  the  problems  of  object  detection, 
it  is  recommended  that  further  work  be  devoted  to  the  following  specific 
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task  s: 


1.  Compute  the  contrast  attenuation  coefficients  for  additional  types  of 
surfaces. 

2.  Determine,  for  a  number  of  selected  cases  of  object-background 
combinations,  the  contrasts  which  one  would  observe  at  the  top  of  the 
atmosphere. 

3.  Introduce  additional  types  of  atmospheric  aerosol  models  and 
determine  their  effects  on  apparent  contrasts  as  seen  from  the  top  of  the 
atmosphere. 

4.  Extend  the  measurements  to  shorter  and  longer  wavelengths. 

5.  Extend  the  measurements  to  more  angles  of  incidence  of  the  radiation 
and  to  the  case  of  diffuse  incident  radiation. 

6.  Extend  the  measurements  to  other  samples  and  surface  conditions, 
and  particularly  to  wet  materials. 

Most  of  the  support  work  which  would  be  required  for  these  tasks 
has  already  been  accomplished.  Essentially  all  of  the  computer  programs 
have  been  written  and  actually  run  on  the  machine,  so  the  computation  of 
additional  cases  is  straight  forward.  The  present  instrument  would  require 
only  relatively  minor  modifications  for  the  additional  measurements,  and 
the  personnel  are  now  completely  checked  out  on  its  operation.  An 
extension  of  the  research  without  interruption  is  strongly  recommended. 
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